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Abstract 



We have measured the transport properties of two mesoscopic hybrid loops 
composed of a normal-metal arm and a superconducting arm. The samples 
differed in the transmittance of the normal/superconducting interfaces. While 
the low transmittance sample showed monotonic behavior in the low tem- 
perature resistance, magnetoresistance and differential resistance, the high 
transmittance sample showed reentrant behavior in all three measurements. 
This reentrant behavior is due to coherent Andreev reflection at the nor- 
mal/superconducting interfaces. We compare the reentrance effect for the 
three different measurements and discuss the results based on the theory of 

quasiclassical Green's functions. 
74.50.+r, 74.80.Fp, 73.23.-b 
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I. INTRODUCTION 



The resistance of a normal metal (N) in contact with a superconductor (S) is modified in 
the vicinity of the N-S interface, a phenomenon well known as the superconducting proxim- 
ity effect [0] . The microscopic mechanism of the proximity effect is Andreev reflection at the 
N-S interface: An electron in the normal metal with energy less than the gap of the super- 
conductor is Andreev-reflected as a hole, with the concurrent generation of a Cooper pair in 
the superconductor . The consequence of this mechanism is the existence of a supercon- 
ducting correlation in the normal metal. Recently, using a theory based on quasiclassical 
nonequilibrium Green's functions, it was shown that the superconducting correlation is ex- 
pected to decay over a length ^(e) = y^D/v/e, where e is the energy of the electron and Dn 
the diffusion constant of electrons in the normal metal The surprising result was 

predicted that the resistance of the normal metal returns to its normal state value at zero 
temperature and energy, the so-called reentrance effect |3|-||]. 

The physical manifestation of this reentrant behavior can be seen in the transport prop- 
erties which are described by an effective diffusion constant D(e,x), a quantity dependent 
on the energy of the electron e and the position x |6|-§. D(e,x) coincides with its normal 
state value at zero energy, increases and reaches a maximum at an intermediate energy of 
the order of the Thouless energy E c = hD^/L 2 , and coincides with its normal state value 
again at higher energies (L here is the length of the normal metal). As a function of tem- 
perature, the resistance of a diffusive normal metal adjacent to a superconductor shows a 
minimum at a temperature T of the order of E c /kB (where ks is the Boltzmann constant), 
and regains its normal state value as T — > 0. A similar minimum is expected in the differ- 
ential resistance of the normal metal as a function of dc voltage V. If the diffusive normal 
metal is connected to two superconductors with different phases, the resistance is expected 
to oscillate as a function of the phase difference between the superconductors, which can be 
modulated by the application of a magnetic field. The amplitude of the magnetoresistance 
oscillations shows a maximum when T is of the order of E c / k& and vanishes again as T — ► 
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(neglecting electron-electron interaction in the normal metal) All this is strictly valid 

only for high interface transmittances; if the interface transmittance is low, the probability 
of Andreev reflection is correspondingly reduced, and the reentrant behavior in the transport 
properties is shifted to lower energy scales, and may disappear entirely [|]|§. 

Several groups have reported observing this reentrance effect in normal metal- 
superconductor (N-S) or semiconductor-superconductor (Sm-S) structures. In N-S struc- 
tures, Courtois et al. |[0| reported observing magnetoresistance oscillations in a normal Cu 
loop with two superconducting Al islands on either side. No reentrant behavior was ob- 
served, however, possibly because it was masked by the Josephson coupling between the 
superconducting islands. Charlat et al. |L1| observed the reentrance effect as a function of 
both temperature and voltage in a Cu loop adjacent to a small superconducting Al island. 
In a different geometry, Petrashov et al. [12| observed the reentrance effect in the ampli- 



tude of magnetoresistance oscillations in an Andreev interferometer. In Sm-S structures, 



den Hartog et al. HT3| , |H| have reported observing reentrant behavior in a geometry where 
a diffusive two dimensional electron gas is coupled to a superconductor to form a loop. In 
these experiments, den Hartog et al. observed reentrant behavior not only in the resistance 
as a function of dc voltage, but also in the amplitude of magnetoresistance oscillations as 
a function of dc voltage. Toyoda et al. [15|] observed reentrance in the magnetoresistance 



oscillations of a two dimensional electron gas connected to a superconducting loop. While 
the results of these experiments qualitatively agree with the theory, quantitative comparison 
with theory is still not satisfactory, especially at high energies or temperatures close to T c , 
where quantitative predictions are difficult to obtain. 

In this paper, we report detailed transport measurements on two N-S devices as a func- 
tion of temperature T, magnetic field H and dc voltage bias V. Both devices are in the 
form of square loops, with one arm of the loop fabricated from a normal metal (Ag or Au) 
and the remaining three arms from a superconductor (Al). The primary difference between 
the two samples is in the transmittance of the N-S interfaces: One sample (sample A) has 
low interface transmittances (Rb/RN > 1, where Rb is the interface resistance and Rn is 



resistance of the normal metal), while the other sample (sample B) has high interface trans- 
mittances {R^/Rjq <C 1). Both samples show a strong temperature dependent resistance 
R(T), large oscillations in the magnetoresistance R(H), and a differential resistance dV/dl 
which is a function of V below the critical temperature T c of the superconductor. While 
sample A shows monotonic behavior in all three measurements, sample B shows reentrant 
behavior in R{T), dV/dI(V), and amplitude of oscillations in R{H) as a function of T. The 
temperature and energy scales for this reentrant behavior are in qualitative agreement with 
recent theories on Andreev reflection in mesoscopic N-S devices although detailed 

quantitative agreement is still lacking. 



II. SAMPLE FABRICATION AND MEASUREMENT 

An electron beam micrograph of sample A is shown in Fig. 1(a), and schematics of the 
two samples A and B are shown in Fig. 1(b). The samples were fabricated by a conventional 
multi-level electron-beam lithography process. The normal metal was deposited first. After 
a second level of e-beam lithography, the normal metal surface was cleaned by a dc Ar + 
etch and the superconductor (Al) was deposited without breaking vacuum in order to ensure 
good contacts between N and S. Au was used as the normal metal for sample A, while Ag 
was used for sample B. Control Ag and Al wires were coevaporated with sample B in order 
to calibrate film properties. The relevant film parameters are as follows: Au/Ag thickness ~ 
28 nm, Al thickness ~ 37 nm, wire linewidth ~ 0.1 —0.14 fim, normal metal (Ag) coherence 
length ^at(T) ~ 0.23 /im/v^T, superconducting coherence length £ai{T = 0) = 0.31 /j,m, 
and electron phase coherence length L v = 0.9 |im at T=30 mK flnj. The area of the N-S 
interfaces was approximately 0.15x0.15 /zm 2 . The samples were measured in a dilution fridge 
between 30 mK and 1.5 K using a four-terminal ac resistance bridge, with ac excitations in 
the range of 10-100 nA, small enough to avoid self-heating. The four terminal measurement 
configuration is shown in Fig 1(b). For the dV/dl measurements, the dc current was applied 
through the same leads as the ac current. Aside from the transparency of the interfaces, the 
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major difference between the two samples is the greater length of the normal metal arms 
beyond the loop in sample A in comparison to sample B (see Fig. 1(b)). 



III. EXPERIMENTAL RESULTS 

Figure 2 shows the resistance of both samples as a function of temperature. The first 
difference noticeable between the two samples can be seen near T c . For sample B, there 
is a sharp drop in resistance as the sample is cooled through T c . If the N-S interface 
resistances were negligible, one might expect the superconducting arm to short out the 
normal arm of the loop, resulting in a decrease in resistance corresponding to the normal 
state resistance of the loop alone. This is indeed the change in resistance we observe for 
sample B within our experimental error, based on the measured resistivities of the normal 
metal and superconductor. Thus, the interface in sample B appears to be highly transparent. 
Since the resistance of sample A does not show a sharp drop at T c , but only a gradual and 
small decrease as the temperature is lowered, we conclude that the interface transmittances 
for this sample are small. We can also estimate the barrier resistance from the resistivities 
of the Ag, Au and Al films along with the measured resistance of the samples. Based on 
these measurements, the resistance of each N-S interface is ~ 25 Q for sample A and < 0.5 
Q for sample B. 

The second major difference between the two samples is seen in the low temperature 
behavior. When the temperature is reduced from T = 1.2 K to 30 mK, the resistance 
of sample A decreases monotonically. Sample B, on the other hand, eventually shows an 
increase in the resistance, resulting in a minimum in the resistance at ~ 520 mK. This 



is similar to the behavior observed by Charlat et al. |TTJ in their Cu/Al loops, and was 
attributed by them to the anomalous proximity effect in the Cu loop induced by the Al island. 
The temperature at which the resistance minimum R m in occurs is given approximately by 
the temperature at which <^v(e = ksT) is comparable to length of the relevant normal 
region. For sample B, the normal regions that contribute to the low temperature zero bias 
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resistance are the small normal arms outside the loop, which have lengths of 0.15 and 0.55 
fim respectively, since the loop itself has zero resistance below T c . The resistance at the 
lowest temperature is 5.9 Q which corresponds to the resistance of the normal side branches 
at this temperature as noted above. Due to its longer length, the contribution from the 0.55 
fim arm should dominate the electrical transport. The temperature corresponding to the 
Thouless energy for this arm is T = E c /kB ~ 170 mK, a factor of 3 lower than the measured 
temperature of ~ 520 mK. However, it should be noted that for a normal wire with one 
end connected to a superconducting reservoir and the other end to a normal reservoir, the 
minimum is expected to occur at T ~ 5E c /ks [4,6]. For sample A, the absence of reentrant 
behavior is consistent with the fact the interface transmittances in this sample are small. 
In addition, the normal side arms are long, and would not be expected to show reentrant 
behavior in our temperature range. If one considers the normal region in the immediate 
vicinity of the low transmission interface as a highly disordered conductor with a very low 
diffusion coefficient D, the relevant £jv(T) is very short, and hence the reentrant behavior is 
pushed to much lower energies and temperatures p.pi|]. 



A similar difference between the two samples can be observed by examining the magne- 
toresistance oscillations as a function of temperature. Figure 3(a) shows the magnetoresis- 
tance of sample A at a few temperatures below 1 K. Oscillations of a period corresponding 
to a flux h/2e through the loop are observed which persist up to the critical temperature T c 
of the superconductor, and whose amplitude at the lowest temperatures is much larger than 
e 2 /h (in terms of conductance). Similar oscillations are seen in sample B (Fig. 3(b)). The 
presence of magnetoresistance oscillations points to the existence of a quantum interference 
effect involving the doubly-connected loop. The large amplitude of these oscillations rules 
out the possibility of their being due to a normal metal quantum interference effect such 
as weak localization or conductance fluctuations, whose amplitude is typically ~ e 2 /h, and 
points to a coherent interference phenomenon involving charge carriers in the normal and 
superconducting arms of the loop pl|-|2"3"|] . 

Figure 3(c) shows the amplitude of the magnetoresistance oscillations for the two samples 
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as a function of temperature. The amplitude is determined by calculating the power in the 
Fourier transform in the inverse field range corresponding to the area of the loop, over the 
field range +/- 25 mT for sample A and +/- 20 mT for sample B. While the oscillation 
amplitude in sample A shows a monotonic increase as the temperature is decreased, the 
amplitude of the oscillations for sample B displays reentrant behavior with a maximum at 
a temperature of ~ 200 mK. Since the oscillations arise from interference effects around the 
loop, one might expect that the amplitude of the oscillations would be determined by the 
ratio of £at(c = ksT) to half the length L of the normal arm, which is ~ 1.1 /im. At T ~170 
mK, 2^Ar(e = ksT) = L. This is in good agreement with the temperature at which we observe 
the amplitude maximum. For sample A, no such maximum is observed, even though the 
film parameters for the two samples are similar. This again is a consequence of the low 
N-S interface transparencies in this sample. At higher temperatures, both samples show a 
temperature dependence which is well described by a function of the form exp[— aL/£jv(T)], 
as can be seen in Fig. 3(c). This is in contrast to the results of Courtois et al. ||10| , where 
the magnetoresistance oscillations were seen to decay as a power law in temperature. For 
comparison, we also show the best fit to the power law dependence found in Ref. [ID|, which 



does not describe the data well. This difference may arise from the difference in the geometry 
of the samples in the two experiments. 

The differential resistance dV/dl as a function of V of sample A and sample B also show 
differences which are consistent with the difference in the quality of their interfaces. Figure 
4 shows dV/dl as a function of V for both samples at low temperature and bias. Sample B 
again shows reentrant behavior, with a resistance minimum at a bias voltage of ~ 7.25 [iV . 
As in the temperature dependent resistance, only the two normal side arms are expected to 
contribute at low dc bias. E c for the longer arm is 15 /J,eV, and hence the voltage at the 
resistance minimum is smaller than expected by approximately factor of two. This should 
be contrasted with the temperature dependence of this sample which was discussed earlier, 
where the temperature at which the minimum in resistance was observed was larger than 
E c /kB by a factor of three. This discrepancy will be discussed later when we attempt to 



compare these data with the quasiclassical Green's function theory. 

In contrast to sample B, sample A shows only a gradual increase in resistance at zero field 
with voltage, consistent with the behavior seen in the temperature dependence. At a finite 
magnetic field of 225 gauss, the curvature of the peak changes, and dV j dl as a function of 
V shows what appears to be reentrant behavior as a function of voltage, similar to that of 
sample B. However, this change in the curvature is not due to reentrance, which can be seen 
by examining the temperature dependence of the sample in a finite magnetic field. Figure 
5(a) shows the temperature dependent resistance of sample A at four different magnetic fields 
corresponding to 0, 1/2, 1 and 3/2 flux quanta h/2e through the area of the loop. Although 
the curves for half-integral flux quanta are different from those for integral flux quanta, no 
reentrant behavior is observed. We believe instead that the change is curvature is similar to 
the zero bias anomaly behavior observed by Kastalsky et al. ||25|| , which was explained by 
van Wees et al. p6| as arising from suppression of coherent multiple Andreev reflections by a 
magnetic field. Figure 5(b) shows similar data for sample B, where curves for both integral 
and half-integral flux quanta show clear reentrant behavior. At half-integral flux quanta, 
there is a small increase in the temperature T min at which the resistance minimum occurs, 
but the curve for zero and one flux quantum are almost the same. This is in contrast to the 
results of Charlat et al. [JTl ] , who saw a monotonic increase in T m « n as the magnetic field 
was increased which they attributed to the field dependence of the electron phase coherence 



length L v . At finite magnetic field, L v is shorter than at zero field [32||. Since L v defines the 
cutoff length for coherent Andreev reflection, L v corresponds to the effective length of the 
sample, and hence the minimum in resistance as a function of temperature would move to 
higher temperatures as a function of magnetic field. This is clearly not seen in our samples. 



IV. DISCUSSION 



s 



A. Quasiclassical Green's function model 



We shall now attempt a quantitative description of the temperature and voltage depen- 
dences of sample B using the quasiclassical Green's function theory. Our analysis is based 
on solving the Usadel equation P7I for the parametrized pair correlation function 9(e, x) in 
the normal metal 9(e, x) is a function of the energy e and position x. Assuming that 

the electron phase coherence length L v is much longer than the length of the normal metal, 
the Usadel equation can be written in the simplified form 

d 2 9{e,x 



dx 2 



+ 2iesin9(e, x) = 



The current is then given by the equation developed in Ref. 
perfect interface): 
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where D(e) is the energy dependent diffusion coefficient in the normal metal which is given 
in terms of 9(e, x) by 
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The conductance is obtained from Eq. (Q) by taking the derivative with respect to the voltage 
V, G(V,T) = dl/dV. At zero bias, the calculation is simplified and the following formula is 
obtained for the resistance R(T): 

de 1 



R(T) = R 



N 



(4) 



/o 2k B Tcosh 2 2j^f \ J L dxsech 2 [Im9(e, x)} 
Equation ([]]) must be solved subject to the appropriate boundary conditions, which are 
usually specified at the N and S reservoirs [^,|]||. At a N reservior, 9(e, x) =0. At a S 



reservoir, 
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for e < A, and 



9(e, x) = i—ln 
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for e > A, where A is the superconducting energy gap. At the N-S interface, 




d8(e,x) 
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dx 



where <jn,s is the conductivity of the interface, S is the cross section of the wire and G b is the 
conductance of the interface. Finally, at a node where two or more normal wires intersect, 
the boundary condition is determined by a Kirchoff-like equation of the form 



where Si denotes the cross section of the branch i joining the node ||. 

Conceptually, at least, determining the resistance of any arbitrary sample is a straight- 
forward matter: one solves the Usadel equation for the pair amplitude 9(e,x) subject to the 
appropriate boundary conditions, then substitutes the result into either Eq. (0) or (f|). Prac- 
tically, however, the Usadel equation needs to be solved numerically. This is not easy, partic- 
ularly for complicated structures such as our Andreev interferometers, and we make certain 
simplifying assumptions to make the calculation tractable. Fig. 6(a) shows a schematic of 
sample B. As we have noted earlier, we assume that the N-S interface resistances are very 
small, so that the loop resistance is zero, as it is shorted by the superconducting arm. The 
measured normal-metal resistance R is then simply the sum of the two N side branches, R\ 
and i?2- To determine the resistance of these structures in the proximity effect regime, we 
need to solve the Usadel equation in the one dimensional wires on either side of the loop. 
Since the electron phase coherence length L v places an upper cutoff to the pair correlation 
in the normal metal, we take the normal reserviors to be at a distance L v from the super- 
conductor. Finally, our calculations show that the effect of the voltage probes on 9(e, x) in 
the side arms of the structure is very small, and hence we ignore the effect of these probes. 

Figure 6(b) shows the final N-S-N geometry that we simulate based on the procedure 
developed in Refs. If we consider the superconductor to be at zero voltage, the 




(5d) 
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energy e of the quasiparticles in each branch % of the structure is related to the voltage 
drop Vi between the corresponding normal reservoir and the superconductor by e, = eV.. 
However, due to the four terminal nature of our measurements, the voltage that is actually 
measured is the voltage V2-V1 at the voltage probes. Although the potential profile in the 
one dimensional normal wire between the normal and superconducting reservoirs is not 



predicted to be linear p8| , the deviations from linearity are small enough that we can relate 
the voltages measured at each probe to the voltage V at the corresponding normal reservoir 
by a linear scaling of the form Vi = V[Lij 'L v ), where Li is the length of the arm from the 
superconductor to the point at which the voltage probe joins the wire. (L v is the effective 
distance to the corresponding normal reservoir.) Furthermore, to calculate the resistance of 
each branch, we need to use Eq. (Q), but with the integral over the length L restricted to 
L{ for each branch, and the normal state resistance Rn corresponding to the normal state 
resistance Rni for each arm. The total temperature dependent resistance of the sample is 
then the sum of the resistances Ri(T) and R 2 (T). 



B. Temperature dependence 

Figure 7(a) shows Ri(T), R 2 (T), and the sum Ri(T) + R 2 (T) calculated in the limit 
A ^> E c for the temperature regime below 1 K. As expected, the major contribution to the 
resistance comes from the longer (0.55 /j,m) side branch. To calculate 6(e,x) we assume A is 
much larger than all the energies integrated in Eq. (f|). For the range of the temperature of 
interest it is good enough to integrate up to 100E C . Although this is larger than the actual 
gap (A = 32E C using BCS theory), changing A in this regime influences the calculation by 
only a small amount in the region of interest. Figure 7(b) shows the result of the calculation 
for the total resistance using two different values of A. In the regime A ^> 100E C , the curve 
recovers its normal state value faster at high temperature than the curve obtained from the 
regime where A is close to 100E C . At temperatures less than 5E C , however, the two regimes 
do not show a significant difference. Therefore, in our simulations we assume A ^> 100E C . 
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Regardless of the value of A assumed in the calculation, however, it is clear that the 
proximity effect theory does not describe well the experimental data in the high temperature 
(T < T c ) regime, as a comparison with Fig. 2 immediately shows. This is because, near T c , 
the contribution from the superconducting transition in the superconducting arm of the loop 
(which we have ignored so far) must be taken into account. For our control pure Al wire, 
which has a width of ~ 0.15 fim, the superconducting transition is fairly sharp, occurring 
within a range of 5-10 mK. The resistance decrease near T c seen in Fig. 2, if it is indeed 
due to the superconducting transition of Al, is much wider. In quasi-one dimensional pure 
superconducting wires, one mechanism that leads to a finite resistance below the nominal 
superconducting transition is nucleation of phase-slip centers. The resistance broadening 
due to these phase-slip centers is given by the Langer-Ambegaokar (LA) form |2^ 



R = —exp 



T\3/ 2 



(6) 



where a is a parameter associated with the attempt frequency for a phase slip event, and /3 
is a parameter related to the energy barrier for a phase slip event. For pure superconducting 
wires, (3 is typically very large (~ 10 6 K) which confines the broadening of the transition 
to a few millikelvin near T c ||29|| . For our N-S structures, since the experimentally measured 
transition is much broader, (3 is expected to be much smaller. 

Figure 8 shows a comparison between the experimental data (triangles) and the theo- 
retical curve taking into account both the proximity effect and the influence of phase slip 
centers. The contribution from the proximity effect (dashed line) was obtained as discussed 
above. The influence of phase slip centers was determined by fitting the difference between 
the experimental data and the proximity effect contribution to the LA equation (j^), using a 
and f3 as fitting parameters. The resulting contribution is shown as the dashed line in Fig. 
8, with the fitting parameters a = 0.17 QK and (3 = 9.93 K. While the contribution due 
to phase slip centers is significant above 0.6 K, the contribution due to the proximity effect 
shows a strong temperature dependence only below 0.4 K. The sum of the two contributions 
is shown as the solid line in the figure, and shows an excellent fit to the experimental data. 
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We should remark that it is not clear that the theory of phase slip centers should be applica- 
ble at all here, although our samples show strong evidence that nonequilibrium phenomena 
are indeed important near the transition []30]]. Nonetheless, our analysis does point to the 



fact that the value of T min clearly depends on the interplay between the contribution due to 
the proximity effect in the normal wire and the decrease in resistance of the superconductor 
near T c , and consequently, the measured T min will not be simply related to the reentrance 
effect. 

C. Voltage dependence 

A similar analysis can be used to calculate the differential resistance dV/dl = R(V) as a 
function of voltage V. To calculate R(V), one notes the temperature kernel in the integrand 
of Eq. (fj) becomes a step function at T = with the discontinuity centered at e = eV, 
i.e., the contribution to the total current / comes only from |e| < eV ||. The conductance 
dl/dV thus contains a 5-function centered at |e| = eV. This results in a simple formula at 
T = §: 

R(V) =R n (t f L s ech 2 [Im9(e, x)} dx) (7) 

\ LJ0 Je=eV 

This formula is applicable to each branch in Fig. 6(b). In order to calculate the total 
resistance R(V), however, one must take into account that the voltage V across the entire 
sample is the sum of the voltages V\ and V2 across each individual branch, subject to 
the condition that the current through both branches is the same. Figure 9 shows the 
experimental data and the calculated curve based on the procedure outlined above. The 
agreement between theory and experiment is clearly not satisfactory. The minimum in 
the experimental resistance occurs at a voltage of ~ 7.25 fj,V, while the theoretical curve 
shows a minimum at ~ 25 fj,V. Using the same rationale as for the temperature dependent 
resistance, we attempt to add the contribution from phase slip centers using the theory of 
LA |2"5fl. However, with the same parameters obtained in Fig. 8, the curve obtained does not 
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agree with the experimental data, primarily due to the strong exponential dependence of the 
theoretical result on the measuring current ||29|| . There are a number of possible reasons for 
this disagreement. First, although we think it is essential that nonequilibrium effects need 
to be taken into account in discussing the voltage dependence, the theory of LA which was 
developed to discuss phase slip centers in pure superconductors may not adequately describe 
nonequilibrium phenomena in N-S devices. Second, our calculations assume perfect N and S 
reservoirs, which are not realized in the experiment. Although this factor might be expected 
to affect both the temperature and the voltage dependence, indications are that the effect 
on the voltage dependence might be more significant pl| . 



D. Magnetoresistance oscillations 

The loop in our Andreev interferometers is essential for the observation of magnetoresis- 
tance oscillations, and hence cannot be ignored in any calculation of the magnetoresistance. 
This complicates the calculation tremendously. Consequently, we have not attempted to 
numerically solve the quasiclassical Green's function equations in the presence of a magnetic 
field. However, the qualitative behavior can be understood by drawing on our experience 
with other quantum interference phenomena in doubly-connected geometries. For the case 
of weak localization in single normal metal rings, for example, the magnetoresistance oscil- 
lates as a function of magnetic field with fundamental period h/2e f32| . The oscillations are 
suppressed exponentially with the phase coherence length L v , exp^—L/L^), where L is the 
perimeter of the loop. In our Andreev interferometers, quantum coherence is maintained 
in the superconducting arms of the loop. In the normal arm, the oscillation amplitude is 
determined by the phase coherence length £jv- Since the normal arm is connected to the su- 
perconductor on both sides, the suppression of the oscillation amplitude might be expected 
to go as ~ exp[— L/(2£n)], where L is now the length of the normal arm. This exponential 
dependence is what we indeed observe in sample A, and also in sample B at higher temper- 
atures. The reentrance effect we see in the amplitude of the magnetoresistance oscillations 
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in sample B is an indication that these oscillations are dependent on the enhancement of 
the diffusion coefficient in the normal arm of the loop. 

V. CONCLUSION 

In conclusion, we have investigated the reentrance effect in two mesoscopic N-S hybrid 
loops with different interface transparencies. The low transmittance sample showed no 
reentrant behavior, consistent with the fact that the relevant energy and temperature scales 
were shifted to values below our measurement range. The high transmittance sample, on 
the other hand, showed reentrant behavior in R(T), dV/dI(V), and amplitude of magne- 
toresistance oscillations, due to the long range coherence of the electron-hole pairs induced 
by Andreev reflection at the N-S interfaces. A quantitative understanding of the experimen- 
tal results cannot be obtained from the quasiclassical Green's function theory of reentrance 
alone. For a more complete quantitative understanding of the properties of such N-S de- 
vices, we believe it is essential to understand the effect of nonequilibrium phenomena on the 
transport properties, particularly at high bias voltages or near the transition temperature 
of the superconductor. 
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FIGURES 

FIG. 1. (a) Scanning electron micrograph of sample A. The additional gate electrode was kept 
grounded and not used in these measurements, (b) Sample schematics for the two samples. The 
dimensions are indicated in fj,m. The leads used to applied ac currents and measure the voltages 
are also shown in the schematics. For the dV/dl measurements, an additional dc current is applied 
through 1+/-. 

FIG. 2. The normalized resistance R/Rn for samples A and B as a function of temperature T. 
i?Ar=67.5 Vt and 10.3 O for samples A and B respectively. 

FIG. 3. (a), (b) are the magnetoresistance curves R{H) for sample A and B respectively. The 
small offset of H in (a) is due to the residual flux trapped in the superconducting magnet. In (a) 
the curves for T=Wl mK, 203 mK, 400 mK and 1.07 K are shifted up by 4 Q, 8 Q, 12 Q, and 16 
Q respectively. In (b) the curves for T = 97 mK, 199 mK, 491 mK and 600 mK are shifted up by 
0.2 Q, 0.4 $7, 0.7 f2 and 0.8 f2 respectively, (c) Normalized amplitude of the Fourier transform of 
(a) and (b) as a function of temperature. The field range is +/- 25 mT and +/- 20 mT and the 
normalization constant is 0.982 Q and 0.019 f2 for samples A and B respectively. The solid lines 
represent fits to the form aexp(— W 1 / 2 ) at higher temperatures, with a = 3.1, 2.7 and b = 3.4, 
2.2 for samples A and B respectively. For comparison, we also show a best fit to a power law of 
the form a/T as used by Courtois et al. [|l(J, with the values a=0.128, 0.238 for samples A and B 
respectively. 

FIG. 4. Normalized dV/dl as a function of dc voltage V at T=30 mK for sample A (solid curve) 
and B (dotted curve). The voltage is obtained by integrating dV/dl vs. Id c - ac and dc currents 
are applied through 1+/- shown in Fig. 1(b). {dV/dI)N is 67.5 Vl and 10.3 Vt for samples A and B 
respectively. The dashed curve shows dV/dl for sample A at a magnetic field of 225 gauss. 
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FIG. 5. R(T) measured for (a) sample A, and (b) sample B at various values of integral 
and half-integral flux quanta h/2e through the loop. Closed symbols, integral flux quanta, open 
symbols, half-integral flux quanta. The temperature dependent curves were obtained from data 
similar to that of Fig. 3. 

FIG. 6. Schematic of the simulated model, (a) Actual sample, (b) Geometry used in the 
calculation. 

FIG. 7. Simulation of the resistance as a function of temperature based on the model in Fig. 
6. (a) Resistance as a function of temperature for R\, R2 and R = R\ + R2 , assuming A ^> E c . 
(b) Total resistance R as a function of T for different A's. A = 200£' c , W 4 E C for the solid curve 
and dashed curve respectively. 

FIG. 8. Theoretical fit to R(T) combining the predictions of the quasiclassical Green's function 
theory with the phase slip model of Langer and Ambegaokar j^9|, as discussed in the text. The 
parameters used in the fit (referring to Eq. (||)) are a = 0.17 f2K, and j3 = 9.93 K, T c = 1.2 K. 
The dotted curve which represents the prediction of the quasiclassical Green's function theory, is 
shifted down by 0.04 f2 and the dashed curve, which represents the predictions of the phase slip 
model of LA is shifted up by 0.45 O for clarity. 

FIG. 9. Comparison of the theoretical calculation (solid curve) based on the theory of quasi- 
classical Green's functions with the measured dV/dI(V) (triangles). The dashed curve shows the 
effect of adding the predictions of the LA theory using the same parameters as in Fig. 8. 
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